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We describe the superconducting proximity effect taking place in a contact between a noncentrosymmetric su- 
perconductor and a diffusive normal/ferromagnetic metal within the quasiclassical theory of superconductivity. 
By solving numerically the Usadel equation with boundary conditions valid for arbitrary interface transparency, 
we show that the analysis of the proximity-modified local density of states in the normal side can be used to 
obtain information about the exotic superconductivity of noncentrosymmetric materials. We point out the sig- 
natures of triplet pairing, the coexistence of triplet and singlet pairing, and particular orbital symmetries of the 
pair potential. Exploiting the directional dependence of the spin polarization pair breaking effect on the triplet 
correlations, we show how the order relation between triplet and singlet gaps can be discriminated and that an 
estimation of the specific gap ratio is possible in some cases. 



PACS numbers: 



I. INTRODUCTION 

Unconventional superconductors are characterized by a 
pairing state breaking not only the f/(l) gauge symmetry. 
The most famous class of such materials is high-temperature 
cuprates in which anisotropic c/-wave pairing is likely real- 
ized.^ It is also known that electrons forming Cooper pairs 
need not to be bound in a spin-singlet state: Spin-triplet states 
with total spin quantum number S = 1 are also possible. 2 
Several materials have been discovered showing strong sig- 
natures of a p-wave spin-triplet pairing state. 3 A full under- 
standing of pairing states in unconventional superconductors 
is still lacking but it is widely accepted that their origin cannot 
be explained within the Bardeen-Cooper-Schrieffer theory of 
conventional s-wave superconductors. 4 The determination of 
the symmetry of the superconducting order parameter is ex- 
tremely important since it is often the first step to understand 
the mechanisms generating the superconducting phase. 

Among unconventional superconductors a special place is 
occupied by noncentrosymmetric superconductors (NCSs), an 
intriguing class of compounds which has recently generated 
much attention due to its unique properties. 5 In these materials 
the superconducting phase develops in a low-symmetry envi- 
ronment with a missing inversion center. This broken symme- 
try generates a Rashba-type 6 spin-orbit coupling (SOC) and 
prevents the usual even/odd classification of Cooper pairs ac- 
cording to orbital parity, allowing a mixed-parity supercon- 
ducting state. Since according to the Pauli principle even 
(odd) orbital parity states are associated with spin-singlet 
(triplet) states, disregarding for the moment the possibility of 
odd-frequency pairing, 7 singlet and triplet Cooper pairs are 
allowed to coexist in NCSs. 

The interest in these materials does not only stem from 
the lacking of a clear understanding of pairing-state symme- 
try and pairing mechanism. Indeed it has been recently rec- 
ognized that NCSs can possibly be in a topologically non- 
trivial phase. This phase is believed to be characterized by 
line nodes, nontrivial topology of Bogoliubov-quasiparticles 



wavefunctions, zero-energy Majorana modes in vortex cores, 
and gapless edge states, all being stable and topologically pro- 
tected against Hamiltonian parameters small variations'"^ 
Moreover they are predicted to be suitable materials to gen- 
erate and host spin polarized supercurrentsf^Sl 

Since the discovery of noncentrosymmetric CePt3Si, ex- 
hibiting a superconducting phase below T c = 0.75K at am- 
bient pressure,^ the family of NCSs has largely grown 
in number. Relevant examples are CeRhSisP CelrShpSl 
CeCoGe 3 P CeIrGe 3 pSI UlrpH Li 2 Pd T Pt 3 _ A Bp] Mo 3 Al 2 cJSl 
and Y2C 3 P^It is worth noting that superconductivity devel- 
oping at oxide interfaces such as LaA10 3 /SrTi0 3 is intrin- 
sically noncentrosymmetric! 2 ^ We will focus on the cerium 
compounds in what follows since within our formalism they 
will admit a common description. They share the same gen- 
erating point group C\ v lacking a mirror plane normal to the 
c axis and have a qualitatively similar T P phase diagram, 
with a supercond ucting dome partially intersecting an antifer- 
romagnetic phase! 20 ! 26 " 29 ' They differ in the amplitudes and 
patterns of the ordered magnetic moments, they have a dif- 
ferent electronic specific heat Sommerfeld coefficient and ef- 
fective masses. CePt 3 Si is the compound with the largest 
Sommer feld c oefficient, a signature of strong electronic cor- 
relationspSEH and it is the only one to become superconduct- 
ing at ambient pressure. A detailed collection of available 
data on these material can be found in Ref. [5] and references 
therein. The determination of superconducting pairing state 
effectively realized is still an open issue. Several experiments 
point toward an unconventional state with the possibility of 
lines nodes! 32 " 36 ! The interpretation of experimental results 
can nevertheless be controversial since strong electronic cor- 
relations and SOC have been show n 30 ! 31 ! to influence the ex- 
perimental findings!^ 7 -! Indeed a triplet pairing state seems to 
be suggested by a critical magnetic field H C 2 larger than the 
condensation energy, e.g. beyond the Pauli-Clogston limitPil 
There is still controversy about this point since a supercon- 
ducting phase developing close to an antiferromagnetic one 
would naturally correspond to singlet pairing state. In this 
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scenario the large critical magnetic field could be explained 
taking into account the SOC.^The standard route to identify 
a triplet pairing state is to exploit the interplay of supercon- 
ductivity and magnetism, namely the magnetic field direction 
dependence of its pair breaking effect on Cooper pairs. While 
the pair breaking exist for any direction of the field for singlet 
superconductors, for triplet superconductors if the field is ori- 
ented parallel to the triplet equal spin pairing direction its ef- 
fect is null. So measuring the spin susceptibility across super- 
conducting transition can help identify triplet pairing states. 
However, the same experimental techniques used to ascertain 
that the prototypical triplet superconductor SrRuC»4 40 has an 
equal spin pairing state with spin oriented in the ab plane 41 
fail when applied to CePt3SiP^While theoretical calculations 
including SOC predict a zero-temperature susceptibility value 
which is half the value of the one in the normal state when 
the field is in a pair-breaking configuration,^ experimentd^Zl 
find a temperature-independent susceptibility for any field di- 
rection. The impossibility to analyze the triplet pairing states 
with the common techniques is a consequence of the strong 
electronic correlations which enlarge the almost temperature- 
independent van V leck term over the Pauli term in the spin 
susceptibility Another complication is that a magnetic 
field shifts the system in the mixed phase which in NCSs is 
predicted to be a helical vortex phase. 43 We believe that these 
complications can be overcome by separating in space super- 
conductivity and magnetism, and that an analysis of pairing 
states realized in NCSs can be put forth by considering hybrid 
contacts with ferromagnets. In particular analyzing the super- 
conducting proximity effect 44 in such hybrids, the supercon- 
ducting correlations propagating in the proximate region will 
be specific of the pairing state in NCSs but will not suffer from 
the complications described above. 

The strategy of using hybrid contacts, junctions, and tunnel- 
ing spectroscopy 45 to probe the pairing state of superconduc- 
tors has been very fruitful. These techniques can be exploited 
to determine the gap amplitude in conventional superconduc- 
tors and even the phase change or the existence of nodes in 
pairing states of unconventional superconductors. 4647 Many 
theoretical studies e xist on tunnel and Josephson junctions 
with NC S s ! 1 ° > 1 2 > 1 4 l 48 n 54 l Here we address the problem of prox- 
imity effect in hybrid contacts of NCSs with normal and fer- 
romagnetic metals. The most direct way to probe the prox- 
imity effect in experiments is to measure the local density of 
states (LDOS) in the nonsuperconducting region and to look 
at how it deviates as a function of position and energy from 
the normal-state LDOS for the influence of superconducting 
correlations. This can be done by placing contacts at differ- 
ent points of the normal region 55 or with a STM tip capable 
of sweeping it! 5 -^ Our aim is to show that this kind of mea- 
surement can be useful to determine both the orbital and spin 
pairing state of NCSs. In particular we point out the signa- 
tures of triplet pairing, triplet and singlet pairing coexistence, 
and particular orbital pair potentials realizations. We show 
how the order relation between triplet and singlet gaps can be 
discriminated and that an estimation of their specific gap ratio 
is possible in some cases. 

The remainder of this paper is organized as follows. In 



Sec. [II] the system is described and the formalism and equa- 
tions employed are introduced. In Sec. [Ill] we report our re- 
sults for the proximity effect both in normal metals and ferro- 
magnets. We compare different types of NCSs with different 
types of gaps. We also consider ferromagnets with varying 
exchange field amplitude and direction. Our conclusions are 



given in Sec. IV 
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FIG. 1: (Color online) In (a) the system under consideration is 
sketched. It is a contact between a NCS and a diffusive ferromag- 
net (DF) separated by an infinitely thin insulator. In such a setup 
the proximity effect features can be probed by looking at how the 
normal LDOS in DF is modified by the presence of the supercon- 
ductor with, e.g., a position and energy resolved STM measurement. 
In the particular case of NCS both even-frequency spin-singlet and 
odd-frequency spin-triplet superconducting correlations are expected 
to propagate in the normal part of the contact; the latter can be long 
ranged depending on the relative orientation of exchange field in DF 
and d- vector in NCS. In (b) the amplitude and phase of the pair po- 
tentials matrices in two different types of NCSs are considered. Nu- 
merical simulations suggest that in cerium family NCSs df poten- 
tials could be stabilized rather than the more commonly assumed sp 
onesP 



II. MODEL AND FORMALISM 

We consider an effective 2D contact between a clean non- 
centrosymmetric superconductor and a diffusive ferromagnet 
(NCS/DF) separated by a thin insulating barrier (see Fig. [T}. 
The contact is at equilibrium at T = 0. The system lies in 
the xy plane with x being the direction orthogonal to the in- 
terface. The NCS side is considered as a reservoir while the 
DF side is assumed to have a finite width L. We describe the 
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contact within the quasiclassical theory of superconductivity, 
a formalism widely used in the analysis of superconducting 
heterostructures (see Ref. 57 for a general introduction and 
Refs.|58]for its application to superconducting hybrids). The 
approximations behind this model have revealed being able 
to grasp the main low-energy experimental features of het- 
erostructures involving conventional and unconventional su- 
perconductors, both for ballistic and diffusive systems. The 
key consideration of the theory is that whenever the Fermi en- 
ergy is much larger than other energy scales in the system, 
only a small error is carried out by considering only physical 
processes taking place close to the Fermi energy. Within this 
formalism, all the physical information is coded in the quasi- 
classical Green's function defined as 



l(r,k F ,e) 



ft J- 



d^ k G(r,k,e), 



(1) 



where G is the position, momentum, and energy dependent 
Fourier-transformed Gor'kov Green's function in the Wigner 
representation and = k 2 /2m. We have chosen units such 
that h = c = 1 and we use . ? . to denote 4x4 matrices in spin 
(g) Nambu space and for 2x2 matrices in spin space. While 
the momentum amplitude dependence has been integrated out, 
the quasiclassical Green's function still depends on Fermi mo- 
mentum direction, i.e., k/r = (cos(p,sin(p,0), where (p is the 
azimuthal angle in the xy plane (see Fig . [TJ . From now on we 
omit the subscript in the momenta since it is clear that in the 
quasiclassical theory they are always fixed on the Fermi sur- 
face. In nonequilibrium conditions a larger number of Green's 
functions together with proper distribution functions are nec- 
essary. This theory, known as Keldysh formalism, is redun- 
dant in our case and the retarded or advanced 4x4 Green's 
functions matrices are sufficient. They have the general struc- 
ture 



zRA 



f 



RA 



(2) 



where g R ' A and f R,A are the normal and anomalous propagator 
blocks, the latter being non zero only when superconducting 
correlations exist. Each block has spin structure 



(R,A RA\ 

% 8 n 
RA RA I 

% Sll J 

and the tilde operation is defined as 

f' A (r,k f ,e)=/^(r,-k F ,-e)*. 
Moreover the matrices are normalized such that 

~«,A ~«A i 

Z -z =1. 



(3) 



(4) 



(5) 



Since retarded and advanced matrices are not independent in 
what follows we will use only the retarded matrix dropping the 
suffix R. Once the latter is known, the LDOS can be calculated 



where No is the density of states in the absence of supercon- 
ductivity. Measuring the LDOS in the nonsuperconducting 
side of superconducting hybrids is the most direct way to 
probe the proximity effect so we will focus on this quantity 
in what follows. The quasiclassical Green's function can be 
calculated by solving the Eilenberger equation, 59 valid in sys- 
tems with arbitrary impurity scattering, and the Usadel equa- 
tion,^ valid in diffusive systems with strong impurity scat- 
tering. Since we are considering a contact between a clean 
superconductor and a diffusive ferromagnet, we have to solve 
different equations in the x ^ regions and their solutions 
have to be connected at x — through a proper boundary con- 
dition. Our choice to consider different impurity scattering 
regimes follows from the consideration that superconductivity 
in noncentrosymmetric materials (and more generally uncon- 
ventional non i— wave super condu ctivity^) is known to occur 
only in rather clean samplesP^^ 

In NCSs the absence of a center of inversion generates an 
electric field which, according to special relativity, in the rest 
frame of electrons is felt as a magnetic field proportional to 
their orbital momenta, resulting in an effective antisymmetric 
Rashba spin-orbit coupling g k = — g k- The particular form of 
this coupling for a noncentrosymmetric material can be deter- 
mined by symmetry considerations. Moreover the broken in- 
version symmetry prevents the classification of Cooper pairs' 
wave function orbital parity which is allowed to be an admix- 
ture of even and odd parts. As a consequence in NCSs an 
admixture of triplet and singlet Cooper pairs can exist. It is 
worth noting that the Pauli principle allows rather than dic- 
tates the existence of a triplet pair potential in NCSs and that 
ideas to identify it are highly desirable. When considering the 
proximity effect in a NCS contact, both singlet and triplet cor- 
relations can manifest in the non superconducting side. Since 
these become isotropic in a diffusive medium, the triplet cor- 
relations have to be odd in frequency in order to respect the 
overall antisymmetry dictated by the Pauli principle. This ex- 
otic type of correlation can decay in a magnetic medium on 
a scale typical of a non magnetic medium, at odds with stan- 
dard even-frequency singlet correlations which are strongly 
suppressed by an exchange fieldPEll n a NCS the gap matrix 
in spin space has the general form 



A k = A s f (k)ia y + A'f (k)d k • aia, , 



(7) 



where A^''' is the singlet (triplet) component amplitude, 
/( s,; )(k) are structure factors, and the well-known d-vector 
dijEH is used to parametrize the triplet pair potentials. The 
usual notation is used for Pauli matrices with a the vec- 
tor of them. It is convenient to normalize Eq. dTll introduc- 



ing the total gap amplitude Ao = maxTr(^A k A k J /2 and the 
triplet/singlet gap amplitudes ratio r = A'/A s E [0,°°), result - 



tlize Eq. d7]l 
Tr(A k A k )/: 



ing in 



A k = Ao ( -4=* «*> + ^ffldk • aia, ) . (Si 



N(e)/N = -Re [Trfe)], 



With the parametrization employed the relative weight of 

(6) 

triplet and singlet gaps amplitudes r can be tuned without af- 
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fecting the total gap amplitude Ao = \/A' s2 + A' 2 and the ex- 
treme cases of pure singlet and triplet superconductors can be 
recovered for r = and r — > °°, respectively. The value as- 
sumed by the ratio r in NCSs is extremely important since 
Andreev bound states can develop at the surface if r > 1 
(A' > A 4 ), depending on the orbital symmetry of pair poten- 
tials. Indeed it has been recognized that NCSs manifest a rich 
and peculiar bound-state structure! 10 " 12 ' So it would be highly 
desirable to estimate experimentally the value of r or at least 
discriminate between the r ^ 1 regimes. Another important 
point is to probe the structure factors and orbital symmetries 
of pair potentials. Several theoretical studies have pointed out 
the possible singlet and tripl et pair p otentials that are sup- 
posed to be realized in NCSsP^^^ Part of our subsequent 
analysis is aimed toward these points. Similar ideas have 
been developed analyzing Raman scattering, tunneling con- 
ductan ce (without p roximity effect) and Josephson effect fea- 
tures! 10 | 12 | 1 ^ Hl51MI it is widely accepted that the spin struc- 
ture of triplet Cooper pair condensate is such that dk | gk since 
this condition maximizes the superconducting critical temper- 
ature.^The form of the spin-orbit coupling vector gk depends 
on crystallographic point-group symmetries of the material!^ 
We focus here on the most common family of NCSs: the 
cerium compounds. Their generating point group is the tetrag- 
onal point group C\ v corresponding to gk « (— A:y,A 5( ,0)!^The 
resulting triplet pair potential has finite S= 1 , m = ± 1 and null 
S = l,m = components. Moreover An and An have equal 
magnitudes resulting in an unitary triplet state. They occupy 
the diagonal elements of the gap matrix while the singlet ele- 
ments are the off -diagonal ones. 

The k dependence of the pair potentials in NCSs is a matter 
of debate. The most common assumption is that they are the 
minimal spherical harmonics consistent with their spin struc- 
ture, i.e., i— wave for singlet and p— wave for triplet pair po- 
tentials. Several studies suggest that this scenario is realized 
m NCSs. 31 This sp case is obtained in our formulation by 
choosing / s (k) = /'(k) = 1. The associated gap matrix A sp is 
depicted in panel (b) of Fig.[T] The singlet off diagonal func- 
tions have opposite sign while the triplet diagonal ones have 
opposite chirality p x ± ip y . Theoretical studies point out that 
higher spherical harmonics pair potentials are possiblef^The 
singlet pair potential can be d x i_ y i— wave and the triplet can 

be /—wave. The gap matrix in the df case A f is depicted 
in panel (b) of Fig. flland is obtained in our formulation by 
choosing / ( (k) = r~v&) = k x —k y . Other possibilities have 
been explored in the literature such as / s (k) = /'(k) = 2k x k y 
which is related to the LaA103/SrTiC<3 heterointerface-^ and 
/ J (k) = 2k x k y , / r (k) = 1 which is not related to any material 
but has been considered for theoretical interest. 53 Part of the 
subsequent analysis is devoted to the comparison of the sp and 
df scenarios, the two most probable pair potential realization 
in NCSs. 

Several studies of NCSs within the quasiclassical theory of 
superconductivity exist in the literature! 13 ! 54 ! 64 ! 70 One often 
employed simplification is to neglect the spin-orbit coupling 
splitting of the Fermi surface since the Fermi energy is always 
much larger than the spin-orbit coupling. It has been pointed 



out how this simplification gives results qualitatively similar 
to the ones obtained without it j 10 l 12 l The application of quasi- 
classical theory to a system with SOC depends on the strength 
of this interaction with respect to Fermi energy: For strong 
SOC one has to calculate two quasiclassical propagators in 
the two bands diagonalizing the SOC while for weak SOC the 
band splitting can be neglected and SOC can be considered 
as a self-energy. 5 While SOC can appear a strong interaction 
with respect to the superconducing gap, it is certainly weak 
with respect to the Fermi energy in the materials under exam- 
ination^ so we will adopt the second strategy. Since we are 
considering a NCS/DF contact where the NCS is in the clean 
limit, in the region x < its asymptotic quasiclassical Green's 
function g$ can be obtained from the Eilenberger equation 



[ep 3 -Ek + Ak,| 5 ] =0, 
where p3 = diag(l, 1,-1,-1) and 



£k 



gk-g 
o 



o 



(9) 



(10) 



After a proper rotation in spin space 70 Eq. [9] can be solved 
analytically. Its solution for generic ratios of triplet and singlet 
gap amplitudes r and for generic structure factors reads 



gs 



+ !-)> 



(11) 



where 



g± 



( N± 
±ie i W± 



Tie- i<f N± ±ie~ i '*A± 



N± 



-A± 

-N ± 



±ie-"*A± ^i e - lv ?N± 
with normal and anomalous elements 

F 



l+r 2 



±ie"*A± 
±ie iv >N± 
-N± J 



(12) 



Ao/± 



^/(l + r 2 )e 2 



where 



/±=f(k)+/(k)r. 



(13) 



(14) 



(15) 



The calculated Green's function does not depend explicitly 
on SOC magnitude but its structure is dictated by its existence. 
This is a well-known result of the application of quasiclassical 
theory to NCSs. 70 However, the absence of SOC is compen- 
sated in our formalism by the parameter r, since the amplitude 
of SOC has been shown to be related to the triplet gap ampli- 
tude in extended Hubbard model studies!^ 6 -! 

The quasiclassical Green's function g for x > in the DF 
can be obtained from the Usadel equation 

D9 x (R|)+i[(e+i5)p 3 + diag[h-CT 1 (h-g)' r ],|] = 0, (16) 
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where D — VpZ/2, x being the relaxation time associated with 
elastic impurity scattering, 8 is an energy associated with in- 
elastic scattering, e.g., quasiparticle damping, "diag" indicates 
a diagonal block matrix, and h is the exchange field vector in 
DF. The characteristic energy scale of the diffusive side is the 
Thouless energy Ejh = D/L 2 . Equation (jT6]t can be solved 
with two boundary conditions at x = and x = L. While at 
x = L the simple Neumann condition 



d x g\ L = 



(17) 



can be applied, the situation at x = 0, where the solution of 
the Usadel equation has to be connected with the asymptotic 
solution of the Eilenberger equation [Eq. (jTTJ], is far more 
complicated. For contacts with conventional superconduc- 
tors this problem has been solved and a boundary condition 
valid for arbitrary interface scattering between the sides of the 
contact has been derived by Nazarov!^ A simplified version 
valid in the tunneling case had been derived by Kuprianov 
and LukichevPS The case of magnetic interface has also been 
treated. 74 In the case of unconventional superconductors a se- 
rious complication arises since the Green's function in the su- 
perconductor has a k (or <p) dependence which is absent in 
the diffusive normal side since the strong impurity scatter- 
ing isotropizes the Green's function. This problem has been 
solved by Tanaka et al. 75 who have generalized the Nazarov 
boundary condition for unconventional superconductors, both 
in a singlet and a triplet state. An infinitely thin insulating bar- 
rier located at x = gives rise to a finite scattering potential 
H5(x) which lowers the transmission probability T as 



r(q>) = 



4 cos 2 (cp) 
4 cos 2 (cp)+ Z 2 ' 



(18) 



where cp is the angle with respect to the normal to the interface 
for a given trajectory (see Fig. [I]), and Z = 2mH/k is the di- 
mensionless barrier strength. The barrier resistance Rg is then 
defined as 



Rb = 



2R 



/-n/2 rf( P r ( ( P) C0S ( ( P) 



(19) 



where Rq = 4-K 2 /k 2 A is the contact Sharvin resistance, A being 
the constriction area. Introducing F = Rg/Rf, with Rp the 
nonsuperconducting side resistance, the boundary condition 
at the interface can be written as 



rL||<A||o = 2([|| , J B(( P )]), 



(20) 



where (. . .) represents a transmission-mediated angular aver- 
age on a half-Fermi surface, 



</(<P)) 



/5/2 c/( P r ( ( P) COS ( ( P) 



(21) 



and B(cp) depends on the asymptotic Green's function in the 
NCS. Defining 



r'(cp) = r(q>)/ h - r(( P )+2 v /i-r(( P ) 



and 

#±(<P) = 2^(9)± &(n-q>)]i 
B((p) can be written as 

-T' (i+HZ l )+T' 2 g\ HZ 1 H + 



fi(<p) = 



v [Ho^r 1 ] +Hi 1 H + -T /2 g \oH: l H + g\ ' 



(23) 



(24) 



(22) 



In Eq. ( |24] > the fraction between matrices is only symbolic and 
should be understood as a matrix product between the inverse 
of the denominator and the numerator. The Usadel equation 
[Eq. ( fT6) i] and the boundary conditions [Eqs. (17) , ( pO) ] form a 
boundary value problem which cannot be solved analytically. 
In the following section we report results obtained by solving 
numerically the boundary value problem with a finite differ- 
ence code that implements the three-stage Lobatto Ilia for- 
mula and an initial guess for the Green's function is adapted 
self-consistently toward the final solution. 



III. RESULTS 

In order to look at the features of the proximity effect in the 
NCS/DF contact, once the quasiclassical Green's function has 
been obtained we calculate the LDOS in the non supercon- 
ducting side with Eq. |6) as a function of energy and position. 
Notice that our LDOS is normalized with the normal state one 
such that when 7V(e) /No = 1 there is no proximity effect. Our 
main target is to find signatures in the proximity effect that 
can give hints on the nature of superconductivity in NCSs- 
namely, if there is triplet superconductivity, what is the ratio 
of triplet and singlet gap amplitudes?-distinguishing between 
sp and df pair potential symmetries. We use Ao as energy unit 
and the superconducting coherence length t, as the length unit. 
The energy is measured from the Fermi energy. The imagi- 
nary energy id in the Usadel equation [Eq. (fTo*)], besides de- 
scribing finite lifetime effects which are always present in real 
systems, has the advantage of stabilizing the code and facili- 
tating the achievement of convergence. It is set to 8 = O.OlAo 
while other parameters will take several values. 

The formalism developed in the previous section allows us 
to explore limiting cases. The asymptotic solution for the 
Green's function in NCSs [Eq. ( fTT) ] describes pure singlet 
and pure triplet superconductors for r = and r — > °°, respec- 
tively. Moreover by neglecting the exchange field in the Us- 
adel equation, the quasiclassical Green's function in a diffu- 
sive normal metal (DN) can be computed. We start exploring 
this case before considering the effect of ferromagnetism on 
the proximity effect. 



A. Proximity effect in NCS/DN 

In Fig. [2] we show the LDOS as a function of energy com- 
paring sp and df NCSs for several r = A' /A* values (r 
0. 0.5, 1 , 2, °° from bottom to top at e = in both panels). The 
LDOS is calculated in the normal side at the border opposite 
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to the interface, i.e., at x — L (see Fig. [TJ, and its width is 
fixed to L = In a typical setup L ~ 10 nm. The interface 
between the NCS and DN is assumed to be in an intermediate 
transparency regime with T = 0.1 and Z = 2. A change in the 
values of these parameters modifies quantitatively the LDOS 
but not its general features, except for L since the proximity 
effect eventually disappears farthest from the interface. In the 
sp case for r — 0, i.e., without a triplet gap, the typical prox- 
imity effect between an s— wave superconductor and a normal 
metal is recovered. Its signature is the opening of a mini- 
gap^l in the LDOS of the normal metal. Its width depends on 
junction parameters such as interface transparency and nor- 
mal layer length.Elln a NCS with < r < 1, e.g., r = 0.5, the 
minigap is still there but its width is smaller with respect to 
the pure s— wave case. Assuming equal amplitude for triplet 
and singlet gaps, i.e., r = 1, the LDOS is almost flattened at 
its normal-state value and a weak proximity effect exists only 
close to the Fermi energy e = 0. This peculiar behavior can 
be regarded as a signature of a transition to a topologically 
nontrivial phase. Indeed the features of the LDOS change 
abruptly for r > 1. The main feature is a well-defined zero- 
energy peak (ZEP) which is associated with the proximity ef- 
fect from triplet superconductors!^ In the df NCS case the 
LDOS is quantitatively less modified by the proximity effect 
since the pair potentials are not fully gapped. The r depen- 
dence of the LDOS in this case is similar to the sp case but 
the energy dependence in the r ^ 1 regimes differs. Indeed a 
dip rather than a minigap exists in the LDOS at low energy for 
r < 1. This is the typical situation in a [100] contact between 
a high T c cuprate and a normal metal while a [110] interface, 
i.e., an effective dxy rather than a d x 2_ y 2 pair potential, would 

not generate any proximity effect. 77 " 79 The dip becomes nar- 
rower increasing r as long as r < 1. Again there is almost 
no proximity effect for r = 1 and a ZEP develops for r > 1. 
This peak is different from the one in the sp case since shoul- 
ders are present at finite energy giving it a bell shape. The 
shoulders are a signature of additional peaks at finite ener- 
gies reflecting the nonmonotonic dispersion of bound states 
in the d f caseP^l This peculiar type of peak is clearly distin- 
guishable from a standard ZEP in tunneling experiments.^ 
The LDOS features described depend quantitatively but not 
qualitatively on the contact parameters; i.e., the width of the 
minigaps and the width and height of the peaks depend on the 
particular values of L, Z, and T, but do not disappear as long as 
proximity effects take place. We have chosen an intermediate 
transparency interface even if this is not the usual experimen- 
tal situation only for the sake of clarity. The main difference 
in considering a tunneling interface is that the peaks and mini- 
gaps appear more narrow, but their features and the possibility 
of their distinction remain unaltered. At the end of this section 
the more proper case of tunneling interface will be considered 
when we report on the penetration length of superconducting 
correlations in the normal side. From the analysis of the prox- 
imity effect in NCS/DN several considerations can be made. 
For both the sp and df cases the regimes r ^ 1 are clearly dis- 
tinguishable by the low-energy LDOS, i.e., minigap for r < 1 
and ZEP for r > 1 in the sp case and dip for r < 1 and ZEP for 
r > 1. The same differences can be exploited to distinguish be- 



tween sp and df NCSs. Even for r > 1 where they both show 
a ZEP in the LDOS, a distinction is possible since only in the 
df case the peak has shoulders at finite energies due to the 
nonmonotonic dispersion of the bound states. The determina- 
tion of the particular value of r rather than the regime r ^ 1 
appears to be more difficult. Indeed it is clear that the features 
of the LDOS in NCS/DN are mostly dominated by the larger 
gap. In other words a NCS with singlet gap larger than the 
triplet one gives rise to a proximity effect very similar to the 
one from a pure singlet superconductor and vice versa. This is 
a general behavior of NCSs not only manifested in the prox- 
imity effect. However, there are quantitative dependencies on 
the particular value of r assumed in both regimes r ^ 1 such as 
the width of the minigap and the height of the ZEP such that 
in principle a fitting procedure can be employed to determine 
its value from the experimental data. 55 A simpler estimation 
is nevertheless possible by looking at the proximity effect in 
NCS/DF. 
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FIG. 2: (Color online) LDOS calculated in the nonsuperconducting 
side of a NCS/DN contact (see Fig. [T} for different NCSs. Several 
r = A' /A s values are plotted in both panels (r = 0,0.5, 1,2,°° from 
bottom to top at e = 0). L = %, T = 0. 1 , Z = 2 are fixed. The LDOS 
is calculated at x = L. 



B. Proximity effect in NCS/DF 

The analysis of the proximity effect in a NCS/DF junction 
can ease the identification of triplet pairing in NCSs, the dis- 
tinction between different types of NCSs, and the estimation 
of the triplet/singlet gap ratio r. Even if the simpler NCS/DN 
analysis can give sufficient information in this sense, loop- 
holes can exist. The point is that observing ZEPs in LDOS is 
not sufficient to claim the existence of a triplet pairing in the 
superconductor since the peak could be generated by uncon- 
ventional singlet superconductivity or could indeed be a sig- 
nature of triplet correlation but generated from spin-sensitive 
processes at the interface or magnetic impurities rather than by 
a real triplet pair potential in the superconductor. The idea is 
to exploit the pair-breaking effect of the exchange field on the 
Cooper pairs. While a spin-polarized environment tends al- 
ways to destroy singlet Cooper pairs, this is not true for triplet 
Cooper pairs where a polarization parallel to the equal spin 
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FIG. 3: (Color online) LDOS calculated in the nonsuperconducting side of a NCS/DF contact (see Fig.[TJ for different NCSs, gaps ratio values, 
exchange field amplitudes and directions in DF. The NCS/DN case (dotted lines) is plotted for comparison. Short-dashed, long-dashed, and 
solid lines are for Ii/Aq = 0.5, 1,5, respectively. L = £, r = 0.1,Z = 2 are fixed. The LDOS is calculated at x — L. 



pairing direction has no effect at allPSEMI Considering that 
the d-vector of triplet superconductivity is built such that the 
equal spin pairing direction is in the plane orthogonal to it, 
this condition reads h_Ldk- If h has some components par- 
allel to dk pair breaking takes place even for triplet Cooper 
pairs. Indeed the standard route to identify triplet pairing is 
to measure the polarization direction dependence of spin sus- 
ceptibility across the superconducting transition. While this 
technique has been proven successful for S12RUO4, the proto- 
typical triplet superconductor, it fails when applied to NCSs 
due to strong electronic correlations. A possible way to over- 
come this difficulty is to extract the superconducting correla- 
tions by the complicated environment of the NCS and analyze 
how they react to a pair-breaking field in a much simpler en- 
vironment such as a ferromagnet, that is, study the proximity 
effect in NCS/DF. In our case the d-vector lies in the xy plane 
so an exchange field out of the plane of the junction h | z is 
not pair breaking while any finite component in the xy plane is 
pair breaking for triplet correlations. For singlet correlations 
the exchange field is pair breaking for any orientation. Fig- 
ure [3] shows the LDOS calculated in the nonsuperconducting 
side of a NCS/DF contact for different NCSs (the first row for 
the sp type and the second for the df type). The two panels 
of a column differ only in the type of NCS with all other pa- 
rameters being equal. Both r ^ 1 regimes are shown for two 
possible orientations of the exchange field h || z and h_Lz. The 
NCS/DN case (dotted lines) is plotted for comparison and in 
each panel short-dashed, long-dashed, and solid lines differ in 
the amplitude of the exchange field, Ii/Aq = 0.5, 1,5, respec- 
tively. L = t„ r = 0. 1, Z = 2 are fixed and the LDOS is calcu- 
lated at* = L. Panels (a), (£>), (e), (/), show how the exchange 
field modifies the LDOS when the NCSs are mostly singlet su- 
perconductors, i.e., A* = 2A ; . In this case the minigap and the 
dip existing in the NCS/DN case are narrowed for low fields, 
while for fields of the order of the superconducting gap an 



enhancement of LDOS is possible and ZEPs in the df case 
or domes in the sp case can develop. These different shapes 
can be used to distinguish between the two types of NCSs. 
The enhancement rather than the suppression of LDOS is in- 
deed found in measurement of contacts between singlet su- 
perconductors and ferromagnets.^ For larger exchange field 
values such as to overcome the Pauli-Clogston limit, the pair- 
breaking effect eventually becomes so strong as to destroy 
completely the proximity effect [see solid curves of panels (b) 
and (f)]. However in a situation of absence of pair breaking 
for the triplet correlations a residual weak proximity effect 
can exist [see solid curves of panels (a) and (e)]. This can be 
exploited as follows. Reporting the NCS/DN case we have 
concluded that a NCS with singlet gap larger than the triplet 
one gives rise to a proximity effect very similar to the one 
from a pure singlet superconductor and vice versa, and that 
a distinction between the extreme cases of pure triplet/singlet 
and mixed cases in the regimes r ^ 1 requires a quantitative 
rather than qualitative analysis. In NCS/DF the situation is 
different. The fact that the exchange field is pair breaking for 
every direction for singlet correlations but not for triplet cor- 
relations gives the possibility to distinguish between a pure 
singlet case (conventional or not) and a real mixed case where 
singlet and triplet gaps exist. Indeed the presence of a sub- 
dominant triplet gap can be ascertained by a residual proxim- 
ity effect in strong ferromagnets with exchange field in a non- 
pair-breaking configuration for the triplet correlations [com- 
pare the solid lines of panels (a), (b), (e), and (f)]. Of course 
in order to observe the pair-breaking effect it is not necessary 
to have a field totally in plane but a finite in-plane component 
is enough. Panels (c), (d), (g), and (h) show how the exchange 
field modifies the LDOS when the NCSs are mostly triplet su- 
perconductors; i.e., A' — 2A S . In this case when the field is 
not in a pair-breaking direction [see panels (c) and (g)], its 
effect on the LDOS is minimal. At low energies there is no 
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difference from the NCS/DN case, while small quantitative 
deviations exist at larger energy. The reason for the insen- 
sitivity of the LDOS to the subdominant component in the 
r > 1 regime lies in the fact that the NCS is in a topologi- 
cally non trivial phase with Andreev bound states at the inter- 
face. These states totally dominate the physics of the contact. 
Once they have been triggered by the condition r > 1, they 
become only weakly dependent on the magnitude of the sub- 
dominant singlet gaps rendering the LDOS almost insensitive 
to an exchange field which is pair breaking only for the singlet 
correlations. 53 In particular the zero-energy bound state with 
momentum orthogonal to the interface does not depend at all 
on the subdominant singlet component and only the bound 
states associated with large-angle trajectories have an appre- 
ciable dependence on it. Moreover the larger contribution to 
the proximity effect comes from almost orthogonal trajecto- 
ries since they are associated with a larger transmission prob- 
ability [see Eq. ([18)], When the exchange field is pair break- 
ing for triplet correlations [see panels (d) and (h)], the peaks 
observed in the NCS/DN case split and lower until eventually 
the proximity effect disappears. The different behavior for 
in-plane and out-of-plane fields can be exploited to ascertain 
whether the peaks observed in the LDOS of NCS/DN are re- 
ally generated by triplet correlations. Indeed only in this case 
a sensitivity on the exchange field direction in NCS/DF should 
be present. 



C. Penetration length of superconducting correlations 

The spatial dependence of the LDOS is worth exploring 
since the penetration depth of superconducting correlations 
in the normal side is also a source of information about the 
pairing state. Figure [4] shows the zero-energy LDOS in the 
nonsuperconducting side of the contact as a function of the 
normalized distance from the interface L/t,. Only the sp sym- 
metry is plotted since at zero energy the df symmetry results 
differ only quantitatively. In each panel F = 10, Z = 10 are 
fixed and a dashed line representing the absence of proximity 
effect is plotted. In panel (a) the case of NCS/DN is shown 
for several r values, r = A'/A v = 0,0.5, 2, °° from bottom to 
top (the first two curves are superposed). It is clear that the 
range of proximity effect in NCS/DN does not depend on the 
triplet/singlet gaps ratio, but when the triplet gap is dominant 
the magnitude of the proximity effect is much more sensitive 
to the distance from the interface. The range of the prox- 
imity effect can have an r dependence in NCS/DF. Panels 
(b) and (c) are for a NCS/DF contact with a mostly singlet 
NCS and with an exchange field out of plane and in plane 
and h/Ao = 5. The cases r = 0,0.4,0.6,0.8,0.9 are plotted 
from bottom to top close to the interface. The oscillation of 
LDOS around the normal-state value and the faster decay with 
respect to the nonmagnetic case typical of singlet supercon- 
ductor/ferromagnet contacts are present for both orientations 
of the field. However when the field is not pair breaking on 
the subdominant triplet correlations an r dependence of su- 
perconducting penetration length clearly appears. The larger 
the triplet component the larger the range of the proximity ef- 
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FIG. 4: (Color online) Zero-energy LDOS calculated in the non su- 
perconducting side of a NCS/DN contact (a) and of a NCS/DF con- 
tact [(b)-(e)] as a function of the normalized distance from the in- 
terface L/t,. Here only the sp symmetry is plotted since at zero en- 
ergy the df symmetry results differ only quantitatively. Everywhere 
r = 10, Z = 10 are fixed and a dashed line representing the absence 
of proximity effect is plotted. In (a) r = A' /A s = 0,0.5, 2, °° are con- 
sidered from bottom to top (the first two are superposed and practi- 
cally indistinguishable). In (b) and (c) the case A' < A s is considered 
for two orientations of the exchange field in the non superconducting 
side and Ii/Aq = 5. The cases r = 0,0.4,0.6,0.8,0.9 are plotted from 
bottom to top close to the interface. In (d) and (e) the case A' > A 1 is 
considered for two orientations of the exchange field in the non su- 
perconducting side and Ii/Aq = 5. The cases r = 1.5, 2, °° are plotted 
from bottom to top close to the interface and everywhere the LDOS 
is calculated at x = L. 



feet. Since the penetration length is easy to measure and the 
quasiclassical theory is able to give a good estimation for it 
from the contact parameters in a conventional superconduc- 
tor/ferromagnet structured the range of the proximity effect 
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could be a precious source of information for NCS supercon- 
ductivity since an unexpectedly long range proximity effect 
would be a strong signature of the existence of a subdomi- 
nant triplet component. This could also be ascertained in the 
NCS's Josephson junction showing a coupling on distances 
larger than expected. When the NCS is mostly a triplet su- 
perconductor this sensitivity disappears as a consequence of 
Andreev bound state formation as explained before. Panels (d) 
and (e) shows this case with an exchange field out of plane and 
in plane and /i/Ao = 5. The curves plotted are for r = 1 .5 , 2, °° 
bottom to top close to the interface. In this case the proximity 
effect is long ranged for out of plane field. Indeed it decays on 
the same scale of a NCS/DN contact [see panel (a)]. In a pair- 
breaking field configuration oscillations around the normal- 
state LDOS appear and the penetration length is reduced as in 
a NCS/DF contact with r < 1 [see panel (c)]. The r depen- 
dence in the r > 1 regime is weak and only small quantitative 
differences emerge close to the interface. One could try to 
overcome this by analyzing the range of correlations at finite 
energies rather than zero energy; however the r dependence 
that emerges is always small compared to the one in the r < 1 
regime. This could be seen by the small deviation induced by 
h in panel (c) and (g) of Fig. [3] From the analysis of the prox- 
imity effect in NCS/DF several considerations can be made. 
The interplay of ferromagnetism and superconductivity can 
give extra information with respect to the NCS/DN case. The 
sp and df types can be distinguished; real equal spin pair- 
ing correlations generated by a triplet gap in NCS can be dis- 
criminated by directional dependence of the pair-breaking ex- 
change field in both energy and space resolved analysis; in 
NCSs with a dominant singlet gap the existence of a subdom- 
inant triplet gap and an estimation of triplet/singlet gap ratios 
are possible. 



IV. CONCLUSIONS 

We have studied the proximity effect in a 2D contact be- 
tween a clean noncentrosymmetric superconductor and a dif- 
fusive metal and ferromagnet within the quasiclassical theory 
of superconductivity. To look at the features of the prox- 
imity effect we have analyzed how the presence of the non- 
centrosymmetric superconductor influences the LDOS in the 
nonsuperconducting side by numerically solving the boundary 
problem determined by the quasiclassical Usadel and Eilen- 
berger equations. We have found unique signatures in the 
proximity effect of the exotic superconductivity expected in 
NCSs and our findings can be easily exploited in tunnel- 
ing spectroscopy experiments, overcoming the difficulties en- 
countered in spin susceptibility measurements. Focusing on 
the most common family of noncentrosymmetric supercon- 
ductors, the cerium compounds, we have found that the ex- 
istence of a dominant triplet gap can be clearly ascertained 
by a zero-energy peak in the LDOS measured in a proximate 
normal metal. We have also shown that the features of the 
LDOS are mostly dominated by the larger gap, a noncen- 
trosymmetric superconductor with singlet gap larger than the 
triplet one giving rise to a proximity effect very similar to the 



one from a pure singlet superconductor and vice versa. How- 
ever, quantitative differences exist and the particular value of 
the triplet/singlet gap specific ratio could be determined by a 
fitting procedure of experimental data. Unique signatures of 
orbital symmetries exist for both triplet and singlet dominat- 
ing gaps. So the proximity effect can distinguish between sp 
and df orbital symmetries. For a dominant singlet gap the sp 
symmetry induces a minigap and the df a dip in the normal 
metal. For a dominant triplet gap both induce a zero-energy 
peak but in the df case it is accompanied by weaker peaks 
at finite energies. The range of superconducting correlations 
penetration in the normal metal does not depend on the or- 
bital symmetries and on the gap amplitude ratio. The analysis 
of proximity effect in a ferromagnet gives additional insights. 
We have shown that it can help to ascertain whether the peaks 
seen in the normal metal case are really a signature of a dom- 
inant triplet gap in the noncentrosymmetric superconductor 
and not a spurious effect generated by magnetic impurities or 
spin-sensitive processes at the interface since a genuine triplet 
zero-energy peak is insensitive to an exchange field oriented 
out of the junction plane while it is destroyed by a sufficiently 
strong in-plane field. In order to change the orientation of the 
field, two separate samples where the exchange field is locked 
to different orientations via, e.g., antiferromagnetic coupling, 
or with different crystallographic orientations can be grown to 
effectively change the orientation of the exchange field. We 
have shown that another signature of triplet dominant gap is 
that the zero-energy peak survives far inside the non super- 
conducting side despite the exchange field when it is in a non- 
pair-breaking configuration. The exchange field can be ex- 
ploited again to distinguish between sp and df types of non- 
centrosymmetric superconductors since for a field of the order 
of the superconducting gap they show different low-energy 
behaviors, sp showing a dome and df a peak in the regime 
where the singlet gap is dominant. We have shown that in 
this regime the presence of a subdominant triplet component 
can be ascertained by the existence of a residual proximity ef- 
fect in strong ferromagnets with exchange field in a non-pair- 
breaking configuration for the triplet correlations. This effect 
can be easier to observe by looking at the penetration length 
of superconducting correlations in the ferromagnet since the 
larger the subdominant triplet component the larger the range 
of the proximity effect. Since this quantity is easy to measure 
and the quasiclassical theory is able to give a good estimation 
for it from the contact parameters in a conventional supercon- 
ductor/ferromagnet structure, the range of the proximity effect 
could be a precious source of information for NCS supercon- 
ductivity since an unexpectedly long range proximity effect 
would be a strong signature of the existence of a subdominant 
triplet component. The specific value of r can be deduced by 
the enhancement in the penetration length in this case. The 
identification of a subdominant singlet component in a mostly 
triplet noncentrosymmetric superconductor appears more sub- 
tle instead since the the Andreev bound states existing in this 
regime are only weakly dependent on the magnitude of the 
subdominant singlet gaps rendering the LDOS almost insen- 
sitive to an exchange field pair breaking only for the singlet 
correlations. 
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The influence of interface scattering on the triplet pair po- 
tential of NCSs deserves a comment since it is expected to be 
pair breaking for anisotropic superconducting order parame- 
ters!^! This effect cannot be easily incorporated in our model 
since we are solving different equations for each side of the 
contact and they are connected by a boundary condition con- 
taining an integration. However, it has been shown 13 that tak- 
ing into account this effect the surface Andreev bound states 
do not disappear but merely change dispersion for large-angle 
trajectories. Since the contribution to the proximity effect of a 
trajectory is inversely proportional to its angle, that would in- 
troduce only small quantitative changes and does not affect 
our conclusions. With the ongoing activity of characteriz- 



ing exotic superconducting materials where triplet and mixed 
pairing are believed to be present we believe our results on 
proximity effect may serve as a useful tool to experimentally 
identify superconducting properties of NCSs. 
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